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CO ' We consider how the presence of strong exciton-phonon coupUng can modify the nature of the 

transition to a condensed phase in systems of strongly coupled light and organic molecules, such 
■ as anthracene and J-aggregates. Considering a model of two-level systems coupled to light, with 

(N ■ relatively strong coupling to phonons, we show that the phonon coupling may drive the condensation 

, transition first order, and that re-entrant behavior vs temperature may arise. We also explain 

how phonon replicas of the polariton spectrum can exist below the chemical potential without 
' macroscopic occupancy. 
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Mierocavity polaritons (superpositions of photons and 
excitons) are light mass bosonic quasiparticles that have 
been observed to form a Bose-Einstein condensate, and 
the non-thermalized cousin, a polariton laser [J, Q . The 
onset of coherence ranges from around 20K in CdTe Q 
to room temperature for materials such as GaN and ZnO 
with large exciton dipole moments : as well as organ- 
ics such as anthracene 0, Q or J-aggregates [§] . Organic 
molecules have localized (Frenkel) excitons that naturally 
have strong optical coupling, but additionally have strong 
coupling to localized molecular vibrational modes, which 
introduces a new set of behaviors and possibilities not 
normally present in inorganic systems. These vibronic 
couplings are mostly studied only at the single exciton 
level, but for example have recently been shown to play 
a complex and crucial role in energy transfer in light har- 
vesting complexes l3"12|. Organic polariton condensa- 
tion provides a natural venue to explore many-particle 
effects arising from simultaneous strong exciton-phonon 
and exciton-photon coupling. In addition, since polari- 
ton condensates and polariton lasers can be seen as an 
exotic variety of ultra-low threshold laser, understand- 
ing polariton lasing with organic molecules may provide 
a pathway to improve the performance of solid state or- 
ganic lasers [l3| . 

The theory of excitons in molecular crystals has been 
extensively developed by Davydov [14| and more recently 
Agranovich [isj . However, modelling organic molecules 
with strong coupling to light is less developed. Work 
to date has focused on using rate equations to calcu- 
late the luminescence spectrum |16| and relaxation pro- 
cesses IrHljas well as the effects of disorder on the spec- 
trum [18|,ll9[. Focusing on low densities, these works have 
generally used a theory of a weakly interacting gas of po- 
laritons, in some cases explicitly derived [2^ from models 
of molecules as saturable absorbers. Such theories well 
describe the very low density regime, but are less accurate 
at higher densities. Instead, one may deal with a model 
of molecules as saturable absorbers (two-level systems) 
coupled to radiation in an optical cavity, by extending 



the Dickc 1211 or Tavis-Cummings [221 model [23|, |24 



which allows one to work at higher densities. 

In this letter, we consider the effects of strong molecu- 
lar exciton-phonon coupling on collective behavior. We 
investigate how the phase transition to a superradiant 
state [25 -27 1 at finite density 2^ 24 1 is modified by the 
coupling to phonons. We find that the phase transition 
can become first order for sufficiently strong phonon cou- 
pling. We also investigate an apparent paradox raised 
by the existence of phonon replicas in the luminescence 
spectrum: Bose-Einstein condensation occurs when the 
chemical potential reaches the bottom of the excitation 
spectrum, at which point a mode becomes macroscopi- 
cally occupied. However, at finite temperature phonon 
replicas exist both above and below the bare polariton 
frequency, and as these are bosonic modes, it appears 
that condensation of these replicas should occur. We find 
that instead, condensation is avoided as the photon spec- 
tral weight vansihes when the chemical potential crosses 
such modes. Nonetheless, at large exciton-phonon cou- 
pling, the polariton mode which does condense contains 
a strong admixture of different phonon replicas. 

The model we study generalizes the Dicke model [2l| . 
which describes two-level systems (molecules) coupled to 
a single photon mode in the mierocavity. To this we 
add a coupling between two-level systems and vibrational 
modes of the molecules. We thus have: 
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+ l}atan + ^^<(an + 4)}- (1) 

Here tp^ is the creation operator for cavity photon, 
ujc = — Mi where ujc is the cavity frequency. Two- 
level systems are described by the Pauli matrices cr^, 
e = e — fjL where e is the optical transition frequency 
of the molecules, the molecule-photon coupling strength 
is denoted g. The chemical potential fi controls the to- 
tal excitation density (excited molecules and photons). 
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which is conserved [2^. Locahzed vibrational modes of 
the molecules are created by the operators djj , these have 
a frequency fl and their coupling to the optical transition 
is given by the Huang-Rhys parameter S. We assume an 
homogeneous system of N two-level systems, and note 
that the characteristic polariton splitting is given by the 
combination gVN which is intensive since g scales as the 
inverse root of the cavity mode volume. This model ne- 
glects finite momentum photon modes, as well as Forster 
transfer of excitations between molecules. The latter is 
unlikely to affect the thermodynamics, as the bandwidth 
due to Forster transfer is small compared to the scale of 
the photon dispersion. Ignoring finite momentum pho- 
ton modes is valid for the high density regime, where 
mean-field approaches (which neglect finite momentum 
fluctuations) become valid [isj . 

In the following we consider the case where the cav- 
ity frequency is detuned above the molecular transition, 
by a value of the order of gVN, i.e. e = cjc — A, A ~ 
gy/N- We focus on this regime for a number of rea- 
sons. Firstly, the clear observation of equilibrated Bose- 
Einstein condensation of polaritons in inorganic materi- 
als 01 required such detuning, so that the lower polariton 
has a strong excitonic fraction, increasing the polariton- 
polariton scattering leading to thermalization [i^ [30| . 
Secondly, the Dicke model can show multiple normal- 
superradiant phase transitions when A > 0; this 
is linked to the existence of "Mott-lobes" in the Jaynes- 
Cummings- Hubbard model [sH . In the Dicke model only 
two such Mott lobes exist. Observation of such behavior 
is not attainable for Wannier excitons, as one cannot eas- 
ily saturate the excitonic occupancy [sl], but may well 
be observable in molecular systems; hence we aim to see 
whether such effects survive coupling to phonons. 

Regarding the vibrational mode frequency we consider 
two cases, fl ~ gVN /2 and fl ^ gVN /20. The first value 
is relevant to anthracene for which gVN ~ SOOmeV 
and the ^ ISOmeV. This is a relatively large vibra- 
tional frequency, due to the stiffness of the fully saturated 
anthracene molecule, and so we also present illustrations 
of the behavior at lower phonon frequency. The value of 
the Huang-Rhys parameter varies significantly between 
different materials; typical values arc around S ~ 1. In 
order to show more clearly the behaviour arising from 
strong coupling, we show below also results that would 
arise with even larger S. 

To calculate the phase diagram of Eq. ([!]) we use 
a semiclassical treatment of the photon field, known 
to be exact 2J] in the thermodynamic limit N — > 
oo, g\/N — > const., i.e. we replace — )- XVN, sub- 
ject to the self-consistency condition UcX = —gy/N{a~). 
Here (cr~) is found by diagonalizing the on-site problem 



+ g^/N{\(T+ + H.c.) + Vta) a 
numerically, truncating the maximum number of vibra- 
tional excitations at Umax 3> S, and thermally populat- 



ing the resulting eigenstates. Anticipating possible first 
order transitions, one must also compare the free energies 
of the normal (A = 0) and condensed (A ^ 0) solutions 
to determine the global minimum free energy. 
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FIG. 1. (Color online). Panel (a) shows the critical value of 
g^N vs (x for various temperatures as indicated. Panel (b) 
shows a map of the order parameter vs temperature and fj, at 
g — 2. Other parameters are A = ojc — e = 4, S = 2, SI = 0.1 

Figure [1] shows phase diagrams, i.e. critical tempera- 
ture or coupling strength as a function of chemical poten- 
tial. As for the regular Dicke model, the behavior at zero 
temperature is reminiscent of normal state Mott lobes, 
separated by a superradiant (condensed) state near 
to fj, = e. At fi — e,T ~ 0, the critical coupling strength 
vanishes. This is analogous to the existence of super- 
fluidity at the boundary between two Mott lobes [33| . 
Coupling to phonons adds several additional features not 
present in the regular Dicke model. Visible in Fig. [Ijb) 
is re-entrant behaviour as a function of temperature, dis- 
cussed below. At larger values of S (see Fig. [2] and later 
discussion) there can be first order transitions into the 
superradiant state. For yet larger S (not shown) a first 
order jump can also occur within the condensed phase. 
(First order transitions have also recently been noted 



within other variants of the Dicke model 

We next discuss these various features in more detail. 
The existence of the re-entrant behavior in Fig. [Ijb) can 
be understood from the effect of phonon sidebands of 
the optical transition meaning that a lower chemical po- 
tential is required to reach the condensed state, i.e. the 
optical transition frequency is reduced due to coupling to 
phonons. For this to happen we require that the tran- 
sition starts from a multiphonon ground state and ends 
in an excited state with fewer vibrations. Hence, such a 
feature is only visible for fc^T > Since the superradi- 
ant bubble vanishes for ksT ^ <?V^, this feature is only 
visible in cases where <^ gVN, Fig. [T{b). 

The existence of a first order normal to superradiant 
transition can be understood by considering the extent 
of polaron formation at large S, i.e. entanglement of the 
phonon and electronic states of the molecule. Such entan- 
glement reduces the overlap between ground and excited 
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states, and so reduces the coupling to light, but also low- 
ers the phonon energy. Such a state is clearly favored 
in the normal state. In the condensed state, it is instead 
favorable to increase the optical polarization by suppress- 
ing the entanglement, and having similar phonon config- 
urations for both electronic states. First order transition 
arises due to a sudden switch between these states. At 
small T the strength of the first order jump is largest at 
points near to, but not exactly at, /.j = e. 

In order to gain a clearer understanding of the first 
order transition, we introduce a variational ansatz which 
captures the behavior at large S. This ansatz can be 
framed as a mean-field treatment of the phonons within 
a coherent state path integral representation. To al- 
low for possibile entanglement between phonons and the 
electronic state, we make a variational polaron trans- 



(a) Exact diagonalization 



(b) Ansatz 



361 on the Hamiltonian H ^ H' 



form [35 

where K = ^y/Sa^' {a) — a). This transformation is exact, 
but if followed by a mean-field approximation d — >■ a, then 
rj, a, and A become variational parameters. This phonon 
mean-field ansatz is valid if the vibrational states are ap- 
proximately coherent states, i.e. when the typical phonon 
occupations, which arc S, are ^ 1. 

The free energy of this phonon mean-field theory can 
be written as: 
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N 
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where C = \/ + (9^)'^ is written in terms of the effective 
molecular transition frequency S and phonon-suppressed 
optical coupling g, given by: 



S = 
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= gVNe'^''''/^. (3) 



Minimizing this free energy with respect to the vari- 
ational parameters A, a and rj we obtain the gap 
equation cZicA ~ g^Atanh(/3C)/2C- Defining k = 
g'^X'^ta.nh{/3Q/[(^ - (J tanh(/3C))^] wc may write the 
equations for a, rj as: 



SVSk tanh(/3C) 

2CK + n '' 
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V 



2CK + n 



(4) 



The equation for rj, describing the extent of polaron state 
formation is instructive. At small A or large 57, r/ — > 1, 
and one has fully developed polarons (vibrational state 
fully entangled with the electronic state, as in the absence 
of any applied field). If, on the other hand, A is large 
then 77 — > and the polaron formation is suppressed; 



i.e. strong driving suppresses polaron formation |35l 
However, an additional level of self-consistency appears 
in the current problem, not present for normal variational 
polaron approaches. The photon field A depends on the 
polarization of the molecules, and the effective coupling 
strength g. Therefore, when the bare coupling g is small, 
the photon field is small, and polarons arc well developed. 
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FIG. 2. Comparison between critical gvN from exact on-site 
diagonalization (left) and phonon mean-field ansatz (right) for 
S = 6, n = 1, A = 4, temperatures as indicated. The color 
scale shows A — {ip)/^/N at the phase boundary. Light (yel- 
low) colors imply second order transition whereas dark (blue) 
first order transition. NB: around resonance the transition is 
always first order, though sometimes weakly 



further suppressing the effective coupling g. At larger 
couplings, polaron formation is suppressed, producing a 
stronger coupling. At zero temperature, this leads to 
a jump within the condensed phase, between a weakly 
and strongly polarized phase. Within the variational ap- 
proach, such a jump occurs near t = ji \i S > 27/8. At 
non-zero temperature the same effect leads to the first 
order normal to condensed transitions, as a competition 
between states with different values of ry. 

In the limit 5 = 0, the phonons are irrelevant, and the 
both phonon mean-field theory and the numerical diago- 
nalization reduce to the Dicke model. Figure [2] shows how 
the two approaches match in the large S limit. In this 
limit, one expects the phonon state can be approximately 
described as a coherent state. However, for smaller S', the 
phonon mean-field theory predicts more strongly first or- 
der transitions than the numerical approach. Numerics 
shows that the first order jump reduces and the transition 
becomes second order as S* — >■ 0. 

When the phase transition is second order, the mecha- 
nism of the phase transition is as for Bose-Einstein con- 
densation: macroscopic occupation of a bosonic mode oc- 
curs when the chemical potential crosses this mode. How- 
ever, this presents an apparent paradox when for a sys- 
tem with phonon replicas. The luminescence spectrum 
shows replicas both above the bare polariton mode, as- 
sociated with creating additional phonons, and at T > 0, 
also below the bare polariton, associated with destroy- 
ing existing thermally populated phonons. Thus at any 
non-zero temperature, an infinite number of (increasingly 
weak) bosonic modes exist below the bare polariton, and 
these appear susceptible to become macroscopically oc- 
cupied at low chemical potential. This clearly does not 
occur in the phase diagrams presented so far. 

The resolution of this apparent paradox lies in the 
changing nature of the luminescence spectrum as one 
varies the chemical potential. The luminescence spec- 
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In a similar manner, one may also find phonon com- 
position of the various modes, and in particular, the 
mode which eventually attains a macroscopic occupation. 
To do this, we may rewrite the above Green's function 
as the inverse of a matrix in the space of photon and 
exciton-phonon states, so that one may determine the 
residue (spectral weight) for a given state, and thus find 
the composition of a given normal mode. Introducing 
Ap,q = g^N (6/3^72-/3^0 _ g-/3£72-/39f2-) ^^^/^g this matrix 
has the form: 



FIG. 3. (a) Spectrum S = 2, T = 0.4, A = 4, J7 = 0.1, 3 = 2, 
where p is the total excitation density, (b) phonon spectral 
weight vs phonon change for the critical mode for various T 



trum arises from the hybridization of the photon with 
the various phonon replicas of the molecular transition, 
and can be found from analytic continuation of the pho- 
ton thermal Green's function [s^l, T^{oJn), given by 



E+_(a;„) 



E 



(5) 
(6) 



where is the photon self-energy, Z is the partition 

function, and apq = |(p|(T~|(7)p is the overlap between 
the polaron states corresponding to p excitations in the 
electronic ground state and q excitations in the electronic 
excited state. The poles of the Green's function cor- 
respond to the excitation spectrum visible in lumines- 
cence, thermally occupied; i.e. the luminescence spec- 
trum P{ljj) = —UBiuj)'^ [-KV^iLUn = UJ — ^ + iO)]. 

The luminescence spectrum calculated from these ex- 
pressions is shown in Fig. |3l^a) , along with the chemical 
potential. It is clear there are many places where the 
chemical potential appears to cross the phonon replicas, 
but crucially, all these crossings are avoided because the 
photon spectral weight of the replica vanishes at these 
points. Even in the regular Dicke model such behavior 
occurs: one can reach an inverted state of the two- level 
systems while avoiding condensation if the chemical po- 
tential crosses the exciton energy at a point where the 
normal modes are pure exciton and pure photon [2^ . 
The same behavior occurs repeatedly in the case with 
phonons. Mathematically this behavior arises from col- 
lision between the zeros and poles of the Green's func- 
tions. The zeros oj* of the Green's function are given by 
uj* —fi = e+{q—p)Q. If i+{q—p)Q = then one may note 
that the numerator in the summand in Eq. ^ disappears 
— at such points poles and zeros cancel, and lj* = 1^1, thus 
whenever a pole and zero coincide, the chemical potential 
crosses a pole. Note however the reverse is not guaran- 
teed; the chemical potential may reach a pole without an 
associated zero, at which point condensation occurs — 
this happens at the right edge of Fig. |3Ja). 
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(7) 



The previous photon Green's function corresponds to 
I?oo, while the various elements T^{pq).{pq) correspond to 
the weights of transitions involving particular initial (p) 
and final (q) phonon number states. 

Figure EKb) shows the composition as a function of 
q ~ p, i.e. the number of absorbed phonons of the mode 
which acquires a macroscopic occupation, exactly at the 
critical point, for various temperatures. At low temper- 
ature, as expected no phonon emission is possible (no 
phonons are thermally populated in the ground state), 
and the distribution depends purely on the overlaps ao,p, 
controlled by the Huang-Rhys factor. For the parameters 
used in Fig.[3l^b), if gVN > fc^T > Q, then re-entrant be- 
havior as shown in Fig. [ijb) occurs. This corresponds to 
the transition at ksT > being associated with macro- 
scopic occupation of a mode involving phonon emission, 
as clearly illustrated in the figure. 

In conclusion, we have shown that strong exciton- 
phonon coupling modifies the phase diagram and the 
spectrum of a polariton condensate, in a parameter 
regime that should be accessible for experiments in or- 
ganic microcavitics. 
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